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We review the phenomenology of the rare decay B — * Xa(.'^ll~ . We present the results of a 
detailed phenomenological analysis and discuss the role of the decay in the search for new 
physics at present and future colliders. Moreover, we extensively elaborate on the size of elec- 
tromagnetic logarithms hi{rnl/m\) in view of experimental cuts. We point out the differences 
^ l' in the analyses of BaBar and Belle and give suggestions on how to treat coUinear photons 

D , in the experimental analyses. These recommendations correspond precisely to theoretical 

y-^ , ' prescriptions and can be combined with measurements performed at a Super-B factory. 

T-H ■ 
^ ' 1 Introduction 

o 

S^ \ In recent years, several flavour facilities as well as the experiments at the Fermilab Tevatron 

have accumulated a large set of data and have confirmed the CKM mechanism [1,2] of quark 

r^ . flavour mixing and CP violation with tremendous success. Besides the determination of the 

^D \ parameters of the CKM matrix and the unitarity triangle, one major goal has been the search 

for new phyics (NP) beyond the Standard Model (SM). Among the prime candidates for this 
search are observables related to rare, flavour-changing neutral current (FCNC) decays of 5, 
D, and K mesons. These decays probe the SM directly at the loop level and are therefore, via 

/\ ' virtual effects, sensitive to scales presently not accessible at direct collider experiments. With 

j^ ■ the start of the LHC being within eyespot, the search for new physics via direct production of 

new degrees of freedom will also become feasible in the near future, and the situation calls for 
an interplay between flavour and collider physics [3-5]. 

Among inclusive flavour-changing neutral current (FCNC) processes (for a review see [6,7]), 
the inclusive B — > Xsl^H." decay represents an important test of the SM, complementary to the 
inclusive B — > Xg^ decay. The two most attractive kinematic observables are the dilepton in- 
variant mass spectrum and the forward-backward asymmetry (FBA), see Fig. [H In the so-called 
'perturbative g^-windows', namely in the low-dilepton-mass region 1 GeV < q^ = mj^ < 6 GeV , 
and also in the high-dilepton-mass region with q'^ > 14.4 GeV , theoretical predictions for the 
invariant mass spectrum are dominated by the perturbative contributions, and a theoretical 
precision of order 10% is in principle possible. 

2 Theoretical framework and Phenomenological results 

The computation of observables of rare decays in flavour physics involves two widely separated 
scales Mh ^ Ml, with Mh — 0{Mw,Mz,mt) and Ml ~ ©(mf,), entailing the need for 
resummation of the occurring large logarithms ln(M^/M£). This is most conveniently done in 



> 

o 



<=■ 1 




Figure 1: Left panel: DifFerential branching ratio (BR) as a function of q^ without (solid, black) and with 
(dotted, red) factorizable cc corrections from the KS approach. The vertical (blue) lines with the arrows indicate 
the perturbative windows. Right panel: Forward-backward asymmetry (FBA) as a function of q^. The meaning 

of the lines is the same as before [9]. 



the framework of an effective theory where the top-quark as well as the heavy electroweak gauge 
bosons are integrated out. In this framework the occurring large logarithms can be resummed 
order by order in a^ by means of techniques of the renormalization group-improved perturbation 
theory. The relevant Lagrangian density can be found, e.g., in Ref. [8]. 

The calculations in i^ ^ Xsi'^i~ have achieved a very sophisticated level. The recently 
calculated NNLL QCD contributions [9-18] have significantly improved the sensitivity of the 
inclusive B -^ Xsi^i~ decay in testing extensions of the SM in the sector of flavour dynamics. 
In particular, the value of the dilepton invariant mass Qq for which the differential FBA vanishes 
is one of the most precise predictions in flavour physics with a theoretical uncertainty of order 5%. 
This well corresponds to the expected experimental sensitivity of 4— 6% at the proposed Super- i? 
factories [19-22]. Also non-perturbative corrections scaling with l/rn^, 1/m^, or 1/m^ [23-29] 
have to be taken into account. Moreover, factorizable long-distance contributions away from the 
resonance peaks are important; here using the Kriiger-Sehgal approach [30] avoids the problem 
of double-counting. 

In the high-g'^ region, one encounters the breakdown of the heavy-mass expansion at the 
endpoint; while the partonic contribution vanishes in the end-point, the 1/mf and 1/mf correc- 
tions tend towards a non-zero value. However, for an integrated high-g'^ spectrum an effective 
expansion is found in inverse powers of ml = rrih x (1 — Vsmin) rather than mi, [31,32]. 

Recently, further refinements were presented such as the NLO QED two-loop corrections to 
the Wilson coefficients whose size is of order 2% [16]. Furthermore, it was shown that in the QED 
one- loop corrections to matrix elements large collinear logarithms of the form log(m^/mwton) 
survive integration if only a restricted part of the dilepton mass spectrum is considered. This 
adds another +2% contribution in the low-g^ region for B{B -^ Xsfi'^fj,~) and results in [8] 
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We will elaborate more on the difference between electron and muon channel in section SI In 
Ref. [33] also the results for the high-g'^ region and for the FBA were derived. The result for the 
branching ratio (BR) in the high-g'^ region reads 
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Figure 2: Left panel: /it-dependence of the forward backward asymmetry for the muonic final state. The lines 
are the NNLO QCD result; the dashed line corresponds to /ib = 5 GeV, and the solid lines to p,h = 2.5, 10 GeV. 
The shaded area is the region spanned by the NLO asymmetry for 2.5 GeV < fit < W GeV. Right panel: Solid: 
FBA as a function of the lepton invariant mass. Curve 2: Reversed sign of C? w.r.t. SM. Curves 1,3: Sign of Cio 

reversed in addition to curves SM,2 respectively [35]. 

In this case the relative impact of the colhnear QED logarithm is about —8% (—20%) for muons 
(electrons) and therefore much larger than in the low-g^ region due to the steep decrease of 
the differential decay width at large q'^. The large error in Eq. ([2]) is mainly due to the sizable 
uncertainties in the parameters that enter the 0{l/mf) non-perturbative corrections. As was 
pointed out in Ref. [29] the error can be significantly decreased by normalizing the B -^ Xsi~^i~ 
decay rate to the semileptonic B — > Xutv decay rate with the same q^ cut. For a lower cut of 
ql = 14.4 GeV^ this leads to 
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where s = q^ /rn^. The uncertainties from poorly known 0{\/m^) power corrections are now 
under control; the largest source of error is Vub- The zero of the FBA is found to be at 

(3.50 ± 0.12) GeV^ 
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The error is about 3% but includes parametric and perturbative uncertainties only. However, 
unknown subleading non-perturbative corrections of order 0{asA/mb), which are estimated to 
give an additional uncertainty of order 5%, have to be added in addition. It is often argued that 
especially the small fi dependence at the zero is an accident and should be increased by hand. 
However, by comparing the NLO-QCD with the NNLO-QCD result one can clearly see that 
the fi dependence is a reasonable reflection of the perturbative error, see left panel in Fig. [2j 
Moreover, the zero is stable under change of the b quark mass scheme; the variation is below 
2% when switching from IS to MS or pole scheme. 



3 New physics sensitivities 

By the end of the B factories the fully differential shape of the branching ratio and FBA will 
not be accessible, contrary to their integrals over bins in the low-g^ region, which are usually 
chosen to be q^ £ [1, 3.5] GeV^ and q^ £ [3.5, 6] GeV^ [34]. These quantities will already allow to 
discriminate between different NP scenarios, see right panel in Fig. [2] as well as Refs. [19,35,36]. 
In the SM the integrated FBA over various bins in the low-q^^ region reads [33] 

(-9.09 ± 0.91)% j(3.5,6) _ r (7.80 ± 0.76)% j(i,6) _ / (-1.50 ± 0.90)% 
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where the upper (lower) hne corresponds to the muon (electron) final state. The relative errors in 
the respective bins are considerably smaller than for the entire low-q^ region since the respective 
values in each bin are similar in size and of opposite sign. The first two numbers do not add up 
to the third one for normalization reasons [33] . In their analysis the authors of Ref . [34] consider 
the three linearly independent quantities that can be extracted from the double differential decay 
width {z = cos d) U 

dPT/{dq^dz) = 3/8 [(1 + z^) Hriq^) + 2 z HA{q^) +2(1- z^) ^^(g^)] , (6) 

where 

dV/dq" = Hriq^) + HUq^) , dA^^/dq^ = 3/iHA{q^) , (7) 

integrated over the aforementioned bins in the low-q^ region. They are able to put contraints 
on the Wilson Coefficients Cg and Cio by imposing a negative Cy whose magnitude is taken 
from B{B — > ^^7). Therefore, if the statistics allows, we highly encourage the experimental 
groups to present their results separately for the three linearly independent observables and for 
the two bins in the low-q^ region. Hence the measurements of the branching ratio and the FBA, 
in addition to the B — > Xg'j branching ratio, will allow to fix magnitude and sign of all relevant 
Wilson coefficients in the SM and to put constraints on the parameter space of NP models. 

New physics might also affect the high-scale Wilson Coefficients in such a way that they 
aquire additional phases. In Refs. [8, 33] we give the results for the branching ratio and FBA 
in terms of generic high-scale Wilson Coefficients. These results may serve to constrain the 
parameter space of NP models and are also of interest in other processes [37]. 

4 CoUinear photons 

After including the NLO QED matrix elements, the electron and muon channels receive different 
contributions due to terms involving ln(m^/m|). This is the only source of the difference between 
these two channels. We emphasize that the results we present in Eqs. ([H)-® correspond to the 
process B — > Xgi'^i^ in which QED photons are included in the Xg system and the di-lepton 
invariant mass does not contain any photon, i.e. q'^ = {pi+ +Pi-Y- This would be exactly the 
case in a fully inclusive analysis using the recoil technique (such an analysis would require a 
Super-S machine [19-22]). 

However, as already pointed out in Refs. [8,33], the presence of the logarithm is strictly 
related to the definition of the dilepton invariant mass. If all photons emitted by the final state 
on-shell leptons are included in the definition of q^: {pi+ + Pi-)^ — > {Pi+ + Pe- +^7)^ then the 
electromagnetic logarithm ln{m'^/m'j) would be absent due to the Lee-Kinoshita-Nauenberg 
theorem [38,39], and hence its effect would disappear Only if all these photons are included in 
the Xg system and not in the q^ (i.e. if a perfect separation of leptons and collinear photons is 
achieved) our expressions containing ln{m'^/mf) are directly applicable. We emphasize that by 
collinear photons we denote only photons that are emitted in the vacuum. No photons emitted 
in the external magnetic field of the detector nor from interaction with the detector material 
are considered. 

In the BaBar and Belle experiments the inclusive decay is measured as a sum over exclusive 
states. Moreover, the treatment of the QED radiation is different [40,41]. As a consequence the 
log-enhanced QED corrections are not directly applicable to the present experimental results 
and have to be modified. Let us elaborate on this. 



°'6 is the angle between the positively charged lepton and the B in the c.m.s. of the lepton pair. 
We note here that if the photons emitted by the final Xs system were also included in the q^ definition then 
there would be an additional ln{Tni/ml). 



For what concerns the di-muon final state, both BaBar and Belle do not include hard collinear 
photons in the q^ definition nor in the final Xg system. So the collinear logarithms associated 
to real photon emission are absent. Hence, the theoretical prediction should contain only the 
logarithms appearing in the calculation of virtual photon exchanges plus real unresolved (i.e. 
soft) emission. In Eqs. ([I])-© we added all real and virtual effects. Because of the details of 
the regularization procedure we used in Refs. [8, 33] (dimensional regularization for both soft 
and collinear singularities), a dedicated calculation is required in order to disentangle soft and 
collinear contributions. 

For what concerns the di-electron final state, both experiments include photons in a cone 
whose opening angle 9 is about (35 — 50) mrad around the lepton directions. The current BaBar 
exclusive, and projected inclusive, st^l~ analyses use only the single most energetic photon 
above the minimum energy threshold (see below) lying in the cone around the electron. For 
photons emitted in this cone we have ?n,| < {pi + p-y)^ < A^ ~ 2E'j{l — cos 9), where Ei is the 
energy of the lepton, usually of order mb/2. Using Ei = mb/2 and = 45 mrad, A is found to 
be of order ?tt,^. Thus, as discussed in Refs. [8,33], if the photons within the cone were included 
in the q^ definition and the collinear photons outside the cone in the final Xs system, our results 
would be directly applicable with ln(m^/mg) replaced by ln(r?2^/A^) where A ~ 0{mp). 

However, the experimental situation is even more complicated. At BaBar and at Belle 
the collinear photons outside the cone are not included at all neither in the q^ nor in the Xg 
contribution. At Belle events with energetic photons (> 100 MeV) or less energetic photons 
(20-100 MeV) are effectively vetoed or suppressed, respectively, by requiring energies and the 
momenta of the final state particles (two leptons, a K^*' and up to four pions) to add up to the 
B meson energy and momentum. At BaBar energetic photons outside the cone are effectively 
vetoed by requiring the momenta of the final state particles (two leptons, a K^*' and up to 
three pions) to add up to the B meson momentum. Moreover, the two experiments differ on the 
definition of the di-electron invariant mass squared, q^: BaBar includes the collinear photons 
inside the cone in the q^, while Belle does not. The first condition requires the subtraction 
of the BR for B -^ Xge'^e'j integrated over photons lying outside of the cone: this effect is 
proportional to log(7TT.^/A^) and can be calculated only by means of a numerical integration over 
the phase space. The second one implies that the log(m^/mg) is absent in the BaBar case and 
present in the Belle one. 

Finally we remind that photons which are resolved but whose energy is below an experiment- 
dependent threshold of about 30 (20) MeV in the lab frame at BaBar (Belle), are not observed, 
i.e., are not included in the di- lepton nor in the Xg system, and effectively shift the energy of the 
reconstructed B mesons. However, most likely these events are still considered to be a candidate 
inside the selection window. These events therefore introduce on the measured branching ratio 
a logarithmic dependence on the aforementioned energy threshold. This dependence can be 
computed by means of the soft photon approximation. In this way all the diagrams factorize 
into tree-level times photon emission. Also the phase space factorizes, and the integral will just 
give the soft and soft-collinear singularities. This kernel can be integrated up to the soft cut 
and when added to the virtual terms, the soft and soft-collinear singularities will drop out and 
will be replaced by the logarithm of the cut. 

The conclusion is that the theoretical predictions for the measurements published by both 
experiments are affected by large collinear logarithms whose precise size has to be determined by 
separate (possibly numerical) calculations which take into account the experimental cuts on the 
photon energy and the opening angle of the collinear cone, as well as the boost from the lab frame 
into the rest frame of the B. We urge both collaborations to converge on a unique definition of 
these processes. We recommend to search and include all photons in the final state and to define 
the di-lepton invariant mass without the photon momentum (i.e. g^ = {p£+ +P£-)'^)'- in this way 
the measured rate and asymmetries would correspond precisely to our theoretical predictions and 



could be combined with future measurements performed with the recoil technique. Alternatively, 
the exclusion of all collinear photons from the measurement of the rate would allow a fully 
analytical calculation of the theoretical prediction (except for of the soft photon cut dependence). 
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